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ABSTRACT 

he  present  lecture  consists  of  two  parts.  In  the  first,  a 
comparison  is  made  of  theory  and  experiment  for  simple  electron 
transfer  reactions  in  solution  and  at  electrodes,  reactions  for 
which  no  bonds  are  broken  or  formed.  The  data  are  numerous  and  were 
recently  reviev/ed  in  part  elsewhere^ 

; The  main  part  of  the  lecture  is  concerned  with  electrode  reac- 
tions involving  the  breaking  and  forming  of  chemical  bonds  and 
principally  with  the  hydronium  ion  discharge  reaction,  H30' 

M(e).<4H20  • H-M(ads.  ).  The  treatment  of  this  reaction  has  been  the 
subject  of  some  controversy  (cf.  1965  Faraday  Discussions  on  Proton 
Transfer  and  more  recent  literature).  To  provide  further  insight, 
a more  dynamical  theory  of  this  elementary  step  for  general  poten- 
tial energy  surfaces  is  described.  It  is  then  shown  how,  depending 
on  the  nature  of  the  potential  energy  surface,  the  theory  leads  to 
various  approaches  in  the  literature.  For  example,  depending  on  the 
surface,  the  system  may  or  may  not  for  dynamical  reasons  be  able  to 
reach  sufficiently  easily  the  saddle-point  region  or.  again,  Franck- 
Condon  solvation  effects  may  or  may  not  contribute  to  the  rate. 
There  is  a clear  need  for  the  application  of  recent  theoretical 
methods  of  electronic  structure  to  determine  the  essential  features 
of  the  potential  energy  surface. 
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ELECTRON  TRANSFER  REACTIONS 


Having  recently  reviewed  the  principal  features  of  the  theory 
of  simple  electron  transfer  reactions  (those  without  bond  rupture 
or  bond  format ion)[l],  I shall,  in  the  interests  of  brev i ty,  call 
attention  to  that  article  and  to  several  other  reviews  [2,3],  The 
notable  contributions  of  Levtch  and  coworkers  [3],  part icularl y, on 
the  quantum  aspects  of  these  electron  transfers  are  well  knovyn. 
Recently,  his  ideas  have  appeared  in  related  approaches  for  under- 
standing the  temperature  dependence  of  electron  transfers  at  very 
low  temperatures  [4]  and  the  behavior  of  highly  exothermic  reactions, 
where  the  process  may  occur  in  the  so-called  inverted  region  [5], 
Theory  and  experiment  in  electron  transfer  reactions  have  had 
a very  fortunate  marriage  during  the  past  two  decades  [1], 

Numerous  comparisons,  correlations,  and  predictions  have  been  pos- 
sible. A major  role  in  the  theory  has  been  played  by  the  Franck- 
Condon  principle  [6];  The  electron  transfers  so  rapidly  chat  the 
other  coordinates  do  not  have  time  to  change  their  values  during 
this  i ns tantaneous  act.  Thereby,  it  was  shown  [7],  their  distri- 
bution must  be  reorganized  before  the  electronic  transition,  if  the 
latter  is  to  occur,  as  well  as  after.  The  relation  between  this 
concept  and  the  role  of  two  intersecting  potential  energy  surfaces 
in  facilitating  Che  electron  transfer  has  been  described  ^8], 

These  simple  electron  transfer  reactions  involve  no  breaking 
or  forming  of  chemical  bonds.  Among  other  things  we  consider  in 
this  paper  the  question  of  when  do  analogous  Franck-Condon  con- 
siderations apply  to  proton  transfers. 


HYDROGEN  ION  DISCHARGE  REACTION 

We  consider  a particular  electrode  reaction  in  which  a bond  is 
ruptured  and  one  is  formed  and  which  has  been  the  subject  of  con- 
siderable attention,  the  hydrogen  ion  discharge  reaction  at  the  in- 
terface between  a solution  and  a metal  electrode  M.  Related  re- 
marks apply  to  proton  transfer  reactions  in  solution.  A mechanism 
for  the  reaction  is  given  by  Eqs,  (l)-(3),  and  we  shall  focus 
attention  on  the  first  step  (1): 


HjO"*  4 M(e)  HzO  + H-M 

0) 

HjO-"  4 H-M(e)  - HjO  4 Hj  4 M 

(2) 

2H-M  - Hj  4 M 

(3) 

The  reaction,  step  (1)  for  example,  has  been  treated  from  two 
very  different  points  of  view  [9-14],  One  view,  used  nowadays 
particularly  by  Bockris,  by  Conway,  and  their  coworkers  [10-14],  is 
an  application  of  transition  state  theory  [14],  wherein  one  cal- 
culates or  assumes  some  potential  energy  surface,  decides  on  the 
transition  state  and  calculates  the  rate.  The  transition  state  is 
typically  assumed  to  be  a hypersurface  passing  through  a saddle- 
point  in  the  potential  energy  surface  or  to  lie  at  the  intersection 
of  a potential  energy  surface  for  reactants  with  one  for  products. 

A second  view,  due  to  Levich  and  coworkers  [9,10],  is  to  assume 
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that  the  individual  or  collective  motions  of  the  reactants  and 
those  of  the  products  and  the  solvent  medium  can  be  treated  as  har- 
monic oscillators,  and  then  use  first-order  oerturbation  theory 
(.Fermi's  Golden  Rule)  or  other  methods  to  calculate  the  reaction 
rate. 

The  relation  between  these  two  approaches  can  be  seen  with  the 
aid  of  a suitable  potential  energy  diagram,  which  will  also  serve 
to  formulate  an  approach,  described  later,  which  in  a sense  encom- 
passes both.  Initially,  we  focus  attention  on  two  coordinates,  the 
0-H  distance  and  the  0-Metal  distance,  and  plot  the  potential  energy 
of  the  system  versus  tfiese  two  coordinates.  Treating  for  the 
moment  the  H?0-H-M  as  a 1 inear  three-center  complex,  one  can  intro- 
duce the  usual  skev/ed  axes  diagram  [lb]  for  three-center  reactions, 
here  HjO.  H and  M.  with  a scaled  0-M  distance  between  the  nearest 
metal  atom  and  the  center  of  mass  of  HjO  as  abscissa  and  a scaled 
distance  between  H and  the  center  of  mass  of  HjO  as  ordinate.  The 
scaling  parameters  and  skewing  angle  are  the  usual  functions  of  the 
three  masses  (i.e.,  of  HjO.  H.  andM),  and  the  kinetic  energy 
then  has  the  same  effective  mass  for  both  scaled  coordinates  [15]. 
Such  a plot  is  given  in  Fig.  1 for  the  case  where  the  reaction  is 
nearly  thermoneutral. 


Fig.  1 Potential  energy  contour  plot  (schematic)  for  the  three 
center  reaction  HjO-H-M  at  a fixed  value  of  the  other 
coordinatesq  and  at  a given  metal -sol ut ion  potential  differ- 
ence. The  mass  of  HjO  is  taken  as  concentrated  on  the  0, 
and  the  mass  of  M as  infinite.  The  rotated  axes  are  scaled 
HjO-HM  and  H-M  distances,  again  with  the  HjO  mass  centered 
on  the  0.  Polar  coordinates  r and  c are  also  indicated. 

The  configurations  along  the  dashed  line  form  the  conven- 
tional transition  state. 


Introducing  polar  coordinates  (r,g)  with  the  apex  of  the  tri- 
angle as  origin,  the  r-motion  corresponds  to  an  0-M  motion  and  the 
Q-mot  ion  to  the  transfer  of  the  proton  between  0 and  M.  /\  profile 
of  the  potential  energy  versus  ra  at  fixed  r is  given  in  Fig.  2 
(solid  curve).  In  addition,  the  potential  energy  is  a function  of 
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all  othei^coord i nates  q,  including  the  various  orientations  of  the 
solvent  molecules  that  solvate  the  three-center  system.  Fig.  2 is 
for  a given  q.  As  r decreases,  the  barrier  height  separating  the 
two  wells  usually  decreases,  and  at  small  enough  r the  double  well 
merges  into  a single  well. 


at  fixed  r and  q,  versus  the  protonic  distance  coordinate 
ra.  Dashed  lines;  a two  electronic  state  ("weak  overlap") 
model  of  this  profile  with  one  curve  for  reactants  R and  the 
other  for  products  P. 


The  difference  in  potential  energy  of  the  minima  of  two 
wells  in  Fig.  2 depends  on  the  electrode-solution  potential  differ- 
ence, and,  so,  thereby,  do  the  barrier  in  Fig.  2 and  the  reaction 
rate. 

Fig.  1 depicts  a potential  energy  surface  for  a roughly  thermo- 
I neutral  reaction.  For  a highly  exothermic  reaction  the  surface 

I resembles  instead  that  in  Fig.  3. 


Fig.  3 Legend  as  in  Fig.  1,  but  the  reaction  is  nav  highly  exother- 
mic. The  dashed  line  again  denotes  the  transition  state. 

I 
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We  can  now  sumriarize  the  two  approaches: 

( 1 ) Transition  Stale  Thet)ry 

Here,  one  typically  expands  the  potential  energy  guadrat- 
ically  in  the  vicinity  of  the  saddle-ooint  of  the  potential  energy 
surface  and  uses  the  local  frequencies  to  calculate  the  partition 
function  and  (tunnel  inyj  transmission  coefficient,  for  use  in  the 
transition  slate  expression  for  the  rate  [15], 

In  the  case  of  a highly  exothermic  reaction  the  transition 
state  is  in  the  reactants'  channel,  as  in  Fig.  3.  In  the  case  of 
a highly  endothermic  reaction  the  transition  state  is,  correspond- 
ingly,  in  the  exit  channel.  In  both  cases  the  reaction  coordinate 
In  the  vicinity  of  the  transition  state  is  seen  to  be  the  0-M  dis- 
tance r whereas  In  the  more  nearl/  thermoneutral  case  In  Fig.  1 the 
reaction  coordinate  in  the  vicinity  of  the  transition  state  was  a 
proton  vibrational  coordinate  a. 

Franck-Condon  effects  on  solvation,  if  any,  are  omitted  in  the 
usual  [15]  transition  state  approach.  The  free  energy  difference 
between  the  transition  state  and  the  reactants  vG  is  computed  and 
from  it  the  reaction  .-ate:  Examples  of  the  application  are  given 

in  [10-14], 

(2)  Radiationless  Transition  Theory  (Levich  and  coworkersj 

The  second  theory  [9.I0]  parallels  that  used  for  radiation- 
less transitions  between  two  et  ect  ronic  states.  The  surface  in  Figs. 

1 and  3 for  reaction  (1)  is  an  el ectron ical 1 v-ad iabat ic  one.  One 
can  imagine  instead  surfaces,  corresponding  to  tv/o  electronic  states 
as  in  Fig.  2 (dashed  lines).  In  one  of  which  the  proton  is  localized 
in  the  reactants'  channel,  i.e. , is  close  to  the  0,  and  in  the 
other  of  which  it  is  localized  in  the  products'  channel,  and  so  is 
close  to  the  M.  Each  local ized  state  has  a series  of  protonic  vi- 
brational energy  levels,  (The  presence  of  a perturbation  forms 
from  the  two  local  ized  electronic  states  two  other  electronic  states, 
of  which  the  lower  energy  one  has  a double  well  solid  curve  profile 
of  F I g.  2.  ) 

Initially,  the  system  (at  r -c)  is  in  one  of  these  localized 
el ec t ron ic -prot on  I c quantum  states.  The  reaction  probability  is 
calculated  at  any  given  r using  Fermi's  Golden  Rule  and  the  above 
perturbation  [9,10],  One  takes  into  account  the  vibrational  over- 
lap integrals  of  the  two  localized  protonic  slates  lust  described 
and  also  of  the  solvation  states  associated  with  each  of  those  two 
protonic  states.  The  appropriate  r for  this  calculation  would  be 
determined  by  some  optimizing  procedure. 

We  also  consider  a third  approach,  one  proposed  here, 

(3)  A Unified  Approach 

(a)  General  considerations 

A detailed  theory  for  reaction  (I)  should  take  into 
account  (i)  the  shape  of  the  actual  potential  energy  surface,  (ii) 
an  appropriate  method,  dynamical  where  necessary,  for  defining  a 
"transition  state,"  and  (iii)  a method  for  handling  the  Franck- 
Condon  effect,  if  any,  on  the  coordinates  q. 

A suitable  ab  i n i t io  potential  energy  surface  or,  in  the  ab- 
sence of  '■hat,  a sem  iemp  i r ica  I surface  such  as  BEBO,  serves  to  cap- 
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ture  at  least  the  qualitative  features  of  (i).  Ir  the  case  of  (ii) 

the  analysis  should  include  the  effect  of  any  large  zero-point  or 

higher  0-H  vibrational  energy  En( r)  in  the  reactants'  channel.  If 

Enlt)  exceeds  the  barrier  in  Fig.  2 at  some  r > r , the 

saddle-point 

H will  transfer  before  the  saddle-point  region  is  reached.  We 
sketch  in  the  next  section  a method  for  determining  a typical  r at 
which  the  proton  transfer  occurs  and  for  determining  dynamically  a 
"transition  state". 

Item  (iii)  is  considered  next,  namely,  anv  Franc k-Condon  effect 
for  the  solvation  coordinates  q.  When  the  oroton  moves  from  the 
reactants'  to  the  products'  channel,  either  near  the  saddle-point  or 
at  a larger  r (for  the  vibrational  energy  reason  just  cited),  it 
moves  so  rapidlv  because  o:  its  small  mass  that  neither  r nor  q nor 

their  conjugate  momenta  change  appreciably  during  that  protonic  motion. 
For  this  reason  a Franc k-Condon  principle  apol ies  to  a proton  jump, 
just  as  it  did  in  electron  transfer  reactions,  though  less  rigor- 
ous 1 y now. 

The  charge  distributions  of  the  tv/o  channels  are  quite  dif- 
ferent but  tend  to  approach  each  other  when  the  two  channels  merge, 
e.  g.  . in  Fig.  I at  the  r corresponding  to  the  saddle-point  region. 

One  can  make  an  estimate  of  the  difference  of  the  charge  distribu- 
tions in  two  channels  and.  thereby,  at  the  most  probable  value  of  r 
for  proton  transfer.  A Franck-Condon  reorgan  izat  i on  is  niseded  when 
the  proton  jumps  between  v;el  1 -separated  channels.  One  can  then 
apply  solvent  reorgan  iza t i on  cons  idenat ions  analogous  to  those  used 
earlier  for  electron  transfer  calculations  M-10].  ho  so' vent 
reorganization  of  this  type  is  needed-for  tra-sfer  at  r's  where  the 
channels  merge, for  the  channels'  charge  distributions  are  then  similar. 

The  effect  of  a "jump"  in  Fig.  1 at  small  r is.  incidentally, 
expected  to  be  smaller  for  proton  than  for  electron  transfers,  since 
in  the  former  case  the  jump  is  probably  usually  a very  short  one 
(e. g. . ~0.  5^)  while  in  the  electron  transfer  case  the  electron 
"jump"  between  charge  centers  which  are  of  the  order  of  5 to  10  ^ 
apart  (typical  ionic  diameters  of  complex  ions). 

In  the  case  of  a highly  exothermic  reaction  t^ig.  5)  depending 
on  whether  the  system  passes  through  the  saddle-point  region  or 
whether  the  proton  jumps  at  a larger  r (because  of  a large  En(r),  say) 
theory  (1)  or  theory  (2)  is  the  more  applicable.  Analogous  remarks 
apply,  by  microscopic  reversibility,  for  the  highly  endothermic 
reaction.  If  the  system  does  pass  through  the  saddle-point  region 
in  the  case  of  Fig.  3,  any  details  of  a subsequent  I,  exothermic)  or 
prior  (endothermic)  proton  jump  are  no  longer  relevant  to  the  cal- 
culation of  the  reaction  rate,  and  conventional  transition  state 
theory  should  again  suffice, 

( b)  Dynamical  Considerations 

To  treat  the  dynamical  question  of  finding  the  r where  most  of 
the  proton  transfer  occurs  the  present  unified  theory  would  make  use 
of  dynamics'  formalism  [16-18]  of  semiclass ical  theory,  and  can  use 
physically-based  approximations  to  simplify  the  dynamics.  The  sys- 
tem is  originally  in  a quantum  state  of  reactants  whose  totality  of 
quantum  numbers  is  denoted  by  n and  proceeds  to  a products'  quantum 
state  m.  The  transition  probability  is  the  square  of  an  S-matrix 
element  It  is  computed  [16-18]  using  classical  trajectories  leading 
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to  products,  beginning  in  state  n and  with  all  possible  initial 
vibrational  phases  of  the  notion.  Whenever  tunneling  occurs,  the 
trajectories  become  complex-valued  rather  than  remaining  real-valued 
[16-18],  and  make  [Smn!^  exponentially  small. 

A dynamical  approximation  which  could  be  adopted  (to  replace 
detailed  traiectory  calculations)  is  that  (a)  the  protonic  oscil- 
lation remains  v i bra t i onal 1 v-ad iaba t ic  [19]  as  long  as  the  trajec- 
tory remains  in  the  entrance  channel  ( i. e.  , the  0-H  vibration  remains 
in  state  n) , (b)  when  each  trajectory  (tunneling  when  necessary) 
crosses  from  the  entrance  channel  to  the  exit  one  it  does  so  essen- 
tially at  a fixed  r because  of  the  rapidity  of  the  protonic  motion, 
and  (c)  thereafter  the  protonic  oscillation  again  remains  vibration- 
al 1 y-ad  i aba  t ic  but  now  in  the  exit  channel  (i.e. , the  H-M  vibration 
is  oresent  in  some  final  state  n).  The  classical  mechanical  analog 
of  a quantum  number  n.  it  will  be  recalled,  is  the  action  variable 
J [19],  related  to  it  by  the  Bohr-Somme rf e 1 d quantization  rule.  e. g. , 
J = (n  - l/2)h  for  a vibration.  Constancy  of  n means,  clissically, 
constancy  of  J.  Classical  adiabaticity  is  discussed  in  [19]. 

Each  initial  phase  w^  of  the  proton  motion  will  determine  the 
value  of  r where  the  H will  strike  the  barrier  in  Fig.  2.  Typically, 
the  system  will  penetrate  the  entrance  channel  as  far  as  it  can  or 
until  the  local  ized  protonic  vibrational  energy  Er|(rj  exceeds  the 
barrier  at  some  r in  Fig.  2.  In  the  former  instance  it  may  reach 
the  products'  channel  without  tunnel i ng  or  by  tunnel ing  at  its  last 
contact  with  the  .^-barrier,  for  the  c-barrier  is  thinner  end  smaller 
there. 


DISCUSSION 

The  preceding  discussion  can  be  used  to  compare  approaches  (l) 
and  (2)  with  each  other  and  with  (3).  Some  differences  in  (1)  and 
(2)  are; 

(i)  The  transfer  in  theory  (I),  in  the  case  of  Fig.  1,  occurs 
in  a saddle-point  region,  usually,  at  small  r.  The  transfer  in  theory 
(2),  based  a zeroth  order  oa i r of  harmonic  oscillators  treatment  plus 
a first-order  perturbation  formal  ism,  is  appropriate  instead  at 
sufficiently  large  r.  (Cf  dashed  lines  in  Fig.  2.) 

(ii)  Since  the  transfer  in  theory  (l),  in  the  case  of  Fig.  1, 
occurs  at  an  r v/here  the  wells  and  hence  the  charge  distributions  of 
the  localized  states  have  merged,  no  solvent  reorgan i zat i on  of  the 
Franck-Condon  type  is  needed.  In  Fig.  3,  since  the  reaction  coor- 
dinate in  theory  (1)  is  the  0-M  distance  any  Franck-Condon  reorgan- 
ization is  irrelevant.  In  constrast,  in  theory  (2)  the  proton  jump 
is  always  rate-determining  and  occurs  at  large  r where  there  is  a 
substantial  separation  of  the  wells  and,  thereby,  a substantial 
difference  of  the  charge  distributions.  The  resulting  Franck- 
Condon  solvent  reorganization  in  theory  (2)  can  be  treated  either  in 
a classical  mechanical  description  [l]  or,  via  Franck-Condon  overlap 
integrals  of  the  collective  solvent  vibrational  wave  functions,  quan- 
tum mechanically  [3]. 

(ill)  The  dynamics  of  the  actual  proton  transfer  is  treated 
differently  in  theories  (l)  and  (2),  usually  being  that  of  a class- 
ical or  a tunneling  motion  in  the  vicinity  of  the  saddle-point  in 
theory  (I)  and  involving  overlap  integrals  and  first-order  pertur- 
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bat  ion  theory  in  theory  (?), 

Both  approaches  represent,  in  a sense,  two  limiting  forms  of 
the  unified  theory. 

Involved  in  the  existing  controversy  was  another  approach  [11] 
in  which  free  energy  for  the  reactants  at  a fixed  OM  distance  v;as 
plotted  versus  a proton  coordinate  [11],  The  free  energy  of  the 
products  was  also  plotted.  The  transition  state  was  assumed  to  1 ie 
at  the  intersection  of  the  two  free  energy  curves.  However,  the  g- 
distribution  tacitly  assumed  for  the  reactants'  plot  vias  the  equi- 
1 ibrium  one  appropriate  to  the  reactants'  charges,  v/hile  that  for 
the  products'  plot  was  a different  equilibrium  one  appropriate  to 
the  products  charges.  Potential  energy  plots  should  be  made  instead, 
and  at  the  same  a-d i s t r i but i on  for  both  curves,  as  discussed  else- 
vjhere  for  homogenous  electron  transfers  [20]  ; 
any  Intersection  of  free  energy  olots  and  the  resulting  discussion 
is  thereby  misleading.  The  Franck-Condon  q-reorgan i zat ion  at  large 
r described  in  the  preceding  section  was  thereby  omitted.  Also  not 
included  was  the  energy  needed  for  penetrating  an  intervening  v;ater 
molecule  by  transferring  H‘  at  large  r's  (~3.  5 R)  or,  if  a 
cavity  was  present,  the  energy  of  first  creating  a cavity. 


CONCLUDING  REMARKS 

During  the  course  of  the  above  discussion  we  have  attempted 
to  answer  a question  which  has  frequently  been  raised.  Does  the 
Franck-Condon  principle  plav  the  same  role  in  proton  transfers  that 
it  does  in  electron  transfers?  When  there  is  a jump  between  two 
distinct  wells,  the  jump  is  so  quick  that  the  solvation  appropriate 
to  the  first  well  in  Fig.  2 does  not  have  time  to  adjust  to  the 
proton  jump  Into  the  second  well.  In  that  case,  there  is  a Frarck- 
Condon  reorganization  of  the  solvation  (the  q-coord i na tes ) . A 
method  was  described  for  determining  the  r v/here  the  proton  "jump" 
occurs.  If  it  occurs  at  an  r where  the  wells  have  merged  the 
Franck-Condon  effect  vanishes.  The  Franck-Condon  is  also  irrel- 
evant in  another  case.  v;hich  was  described. 

Several  theories  were  compared,  the  source  of  their  differences 
analyzed,  and  a "unified"  theory  outlined. 
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